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We investigate the ground-state phase diagram for the ferromagnetic Kondo lattice model 
on a triangular lattice by a variational calculation for various spin orderings up to a four-site 
unit cell. We find that a noncoplanar four-sublattice ordering with a finite scalar spin chirality 
emerges at and around 1/4 filling, in addition to the 3/4-filled case, which was predicted to be 
induced by the perfect nesting of the Fermi surface [I. Martin and C. D. Batista: Phys. Rev. 
Lett. 101 (2008) 156402]. The 1/4-filling phase is stable in a wider range of parameters than 
the 3/4- filling one, and includes a large region of gapped insulating state characterized by a 
Chern number. We also compute the Hall conductivity in the chiral-ordered phases. 
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Spin chirality has been an issue of interest as an emer- 
gent composite degree of freedom in a broad range of 
condensed matter physics. NoncoUinear spin config- 
urations give rise to a finite spin chirality, e.g., the vec- 
tor spin chirality 5i x ^2 -I- ^2 x 5^3 -f 5*3 x Si or the 
scalar spin chirality 5*1 • (52 x 6*3), defined on a triangle 
of spins. The latter scalar chirality, which becomes finite 
for a noncoplanar spin configuration, has recently drawn 
considerable attention not only in localized spin systems 
but also in itinerant electron systems as an intrinsic ori- 
gin of the anomalous Hall effect (AHE). When itinerant 
electrons are coupled to such specific spin configuration, 
electrons experience an internal magnetic field according 
to the solid angle spanning three spins through the so- 
called Berry phase, which can result in AHE. This new 
mechanism was first theoretically proposed for a chiral 
order in the ferromagnetic Kondo lattice model on a two- 
dimensional kagome lattice, and developed in a three- 
dimensional frustrated system.^) The relevance has been 
discussed in experiments on metallic pyrochlore oxides. 

Recently, this chirality mechanism was explored for 
triangular lattice systems. Martin and Batista pointed 
out the possibility of AHE in a Kondo lattice model on 
the triangular lattice at 3/4 filling of itinerant electrons, 
which is brought about by a noncoplanar four-sublattice 
spin ordering due to the perfect nesting of the Fermi 
surface.*' This proposal is attractive since similar non- 
coplanar four-sublattice orders are encountered in some 
triangular-lattice spin systems.^' 

Although previous studies have revealed a nontrivial 
relation between chirality ordering and AHE, they rely 
on the assumption that the spin configuration is a given 
texture and not affected by the coupling to itinerant elec- 
trons. In spin-charge coupled systems, however, the spin 
state is determined in a self-consistent manner with the 
itinerant electron state through the interplay between 
spin and charge. Hence, it remains unclear whether such 
chiral ordering is realized under spin-charge interplay. 
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For further understanding of spin-chirality-related phe- 
nomena, it is necessary to carry out a microscopic anal- 
ysis of the most stable spin state and to clarify whether 
the chiral order with AHE remains stable even when con- 
sidering the interplay between spin and charge. 

In this study, we investigate the ground state of the 
ferromagnetic Kondo lattice model on the frustrated tri- 
angular lattice with emphasis on the possibility of chiral 
ordering. We compare the energies for various spin con- 
figurations up to four-sublattice orders, and determine 
the ground-state phase diagram. As a result, we find 
that a noncoplanar four-sublattice spin ordering with 
a finite scalar spin-chirality emerges around 1/4 filling, 
in addition to the 3/4 filling predicted in the previous 
study.*) This new phase includes both metallic and insu- 
lating phases and is stabilized in a wider parameter re- 
gion than the 3/4-filling phase. We also discuss AHE in 
these chiral-ordered phases by calculating the Hall con- 
ductivity, which is quantized in the insulating phases cor- 
responding to the associated Chern number.^-*' 

We consider the ferromagnetic Kondo lattice model on 
the triangular lattice as one of the fundamental models 
for describing the interaction between localized spins and 
conduction electrons under frustration. The Hamiltonian 
is given by 
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where c|Q,(ci_Q,) is a creation (annihilation) operator for 
a conduction electron with spin a on site i, and t is 
the transfer integral, Jh is the Hund's-rule coupling, 
CTq^ = {^afi^'^^ap^'^aj}) ^ vector of Pauli matrices, Si 
is a localized spin on site i, and Jk is the antiferromag- 
netic superexchange interaction between localized spins. 
We consider classical spins for Si with \ Si\ = 1. The sums 
are taken over the nearest-neighbor sites on the tri- 
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Fig. 1. Ordering patterns of localized spins used in the varia- 
tional calculations: (la) ferromagnetic order, (2a) and (2b) two- 
sublatticc orders, (3a)-(3d) three-sublattice orders, and (4a)-(4f) 
four-sublattice orders. See text for details. 

angular lattice. Hereafter, we take t = 1 as an energy 
unit. 

We investigate the ground state of the model given 
by eq. (1) while varying the electron density n ~ 
W^iai^i.a^i,a) IS the total number of sites), Jh, 
and Jk- We compare the ground-state energies for dif- 
ferent ordered states of the localized spins and determine 
the most stable ordering by a variational calculation. In 
the calculation, we consider 13 different types of ordered 
states, up to a four-sublattice unit cell, as shown in Fig. 1. 
Figure l(la) shows a ferromagnetic order. Figures l(2a) 
and l(2b) show two-sublattice orders: (2a) a collinear 
stripe order and (2b) a stripe order with a canting angle 
9. Figures l(3a)-l(3d) show three-sublattice orders: (3a) 
a 120° noncollinear order, (3b) a noncoplanar umbrella- 
type order with angle 9 (canted in the normal direction to 
the coplanar plane from the 120° order), (3c) a coplanar 
order with canting angle 6 for two spins from 120° order, 
and (3d) a 2:l-type order with two parallel spins that 
have angle 6 to the remaining one. Figures l(4a)-l(4f) 
show four-sublattice orders: (4a) an all-out-type order 
which was discussed in Ref. 8, (4b) a two-in two-out-type 
order, (4c) a three-in one-out-type order, (4d) an all-out- 
type order with canting angle 9 for three spins, (4e) a 3:1 
collinear order, and (4f) a coplanar order with a 90° flux- 
type configuration. For the given parameters n, Jh, and 
Jk, we determine which spin state is most energetically 
stable by comparing the ground-state energies for all the 
spin states listed above. For the states (2b), (3b), (3c), 
(3d), and (4d), we optimize the canting angle 6. Note 
that (2b) with 6* = tt, (3b) with 6* = f , (3c) with 9 = 0, 




Fig. 2. (Color online) (a) Grand-canonical potential at T = 0, 
O = CH) — /in, for relevant ordered states measured from that 
for the ferromagnetic order and (b) the electron density n as 
a function of the chemical potential fi. The data are for Jh = 
2.0 and Jj^ = 0. Phase-separated regions, shaded areas in (b), 
are determined from the jumps of n in (b) associated with the 
crossing points in (a). 



(4d) with 9 = cos-i(-l/3), 9 = cos-i(+l/3), and 6l = tt 
are equivalent to (2a), (3a), (3a), (4a), (4c), and (4e), 
respectively. Later, to check the local stability of (4a), 
we will also try a general configuration with all possible 
spin directions within the four-site unit cell. Although an 
incommensurate order might take place for a general fill- 
ing, we focus on uniform q = orders with the magnetic 
unit cells listed above, particularly, on their stability at 
and around commensurate fillings, as we discuss below. 

Figure 2(a) shows a comparision of the grand- 
canonical potential at T = 0, = ("H) — ^n, for vari- 
ous ordered states as a function of the chemical poten- 
tial fj. for Jh = 2.0 and Jk = 0, as an example. The 
data are computed by approximating the integral over 
the (unfolded) first Brillouin zone using the sum over 
grid points of 1600 x 1600. We here plot for relevant 
ordered states at these parameters measured from that 
for the ferromagnetic state in Fig. l(la). (All other states 
not shown have higher f2.) In the low- and high-/i regions, 
the ferromagnetic state (la) gives the lowest $7. However, 
in the intermediate region, —3.3 < /x < 1.6, other states 
have a lower fl. For each range, the state that has the 
lowest gives the ground state. 

To find the ground-state phase diagram as a function 
of n, it is necessary to examine the relationship between 
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Fig. 3. (Color online) Ground-state phase diagram for (a) Jk = 
and (b) Jk = 0.01. The regions (la)-(4f) correspond to the 
ordering patterns in Figs. l(la)-l(4f), respectively. The vertical 
dashed lines at n = 1/4, 1/2, and 3/4 show gapful insulating 
regions. PS indicates a phase-separated region. 



n and fi. In general, a phase transition between different 
magnetic orders is of first order and accompanies a jump 
of n. This is demonstrated in Fig. 2(b), which plots n as 
a function of /i. For example, at /i ~ —3.3, the ground 
state changes from a ferromagnetic state [Fig. l(la)] to a 
four-sublattice all-out-type ordered state [Fig. 1 (4a)] , as 
shown in Fig. 2(a). At the same fi, as shown in Fig. 2(b), 
the density n has different values for these two states, 
i.e., n ~ 0.16 for the former and n ~ 0.23 for the latter, 
which results in a discontinuous change of n at the phase 
transition. Since the system is unstable and cannot have 
a fixed density in the discontinuous regime, the jump of 
n signals a phase separation. Similar analysis identifies 
phase-separated regions, as demonstrated in Fig. 2. 

We perform the comparison of f2 and the identification 
of phase-separated regions while varying the parameters 
n, Jh, and Jk- For Jh <C t, in particular, at Jk = 
0, the difference in f2 among different orders becomes 
very small. In such cases, we determine the ground states 
using the perturbation theory in Jh /t up to fourth order. 

Figure 3(a) shows the result of the phase diagram as 
functions of n and Jh at Jk = 0. In the low- and high- 
density regions, a ferromagnetic metallic phase appears 
[Fig. l(la)]. The phase becomes wider as Jh increases. 
This ferromagnetic phase is stabilized by the double- 



exchange mechanism. On the other hand, a three- 
sublattice 120° order emerges near half-filling n ~ 1/2 
[Fig. l(3a)]. In particular, the system becomes insulating 
at n = 1/2 for Jh > 2. This is understood by considering 
the fact that, at half-filling, the second-order perturba- 
tion in t/ Jh leads to an effective antiferromagnetic inter- 
action between localized spins. Under the frustration in 
the triangular lattice, the antiferromagnetic interaction 
stabilizes the 120° order. In the large- region, a phase 
separation takes place between the ferromagnetic state 
and the 120°-ordered state. 

For smaller Jh, the phase diagram becomes more com- 
plicated. Among various phases, the most interesting 
point is that a four-sublattice all-out order with a finite 
scalar chirality [Fig. l(4a)] is realized near 1/4 filling and 
3/4 filling. The phase at 3/4 filling is the one predicted 
by Martin and Batista as a consequence of the perfect 
nesting of the Fermi surface.^' Our energy comparison 
confirms the prediction. In contrast, the phase around 
1/4 filling is not of nesting origin, since the Fermi sur- 
face is almost circular at this filling for Jh = 0. This is 
a new chiral phase with a four-sublattice order, which 
was not predicted in the previous work. Because of the 
noncoplanar spin order, this chiral phase exhibits AHE 
as in the 3/4 filling phase; we will return to this point 
later. 

As shown in Fig. 3(a), the new phase at n ~ 1/4 
is stabilized in a wider range of parameters than the 
nesting-driven 3/4-filling phase; the former is stabilized 
for < Jh < 4.5, whereas the latter is limited for 
< Jh < 0.51. Moreover, the 1/4-filling phase becomes 
insulating in a wide Jh range of 0.69 < Jh < 4.5. The 
gap opens between the lowest and second bands of the 
four bands under the four-sublattice order (each band 
has a twofold degeneracy in terms of spin). This gap 
appears to robustly stabilize the chiral order; in fact, 
we confirm that the gapped chiral state remains stable 
even when allowing all possible spin configurations while 
changing the relative angle of spins continuously within 
the four-site unit cell. 

When we switch on the superexchangc interaction Jk, 
the chiral phases as well as the threc-sublattice 120°- 
ordered phase become more stable. Energy changes are 
easily calculated from the last term in eq. (1), for ex- 
ample, —Jk (— 3Ji<-/2) per site for the four-sublattice 
all-out (three-sublattice 120°) ordered state. Figure 3(b) 
shows the phase diagram at Jk = 0.01. Indeed, the 120°- 
ordered state covers a wider range of parameters; in par- 
ticular, it extends to the small- Jh region and n ^ 1 re- 
gion. In addition, the four-sublattice chiral phases around 
1/4 and 3/4 fillings are also enlarged. Instead, the ferro- 
magnetic regions are reduced and surrounded by phase 
separations up to the 120°-ordered state. Thus, Jk sta- 
bilizes the chiral phases as well as the 120°-ordered state. 

Let us discuss AHE for the model eq. (1), while fo- 
cusing on the chiral phases near 1/4 and 3/4 fillings. We 
calculate the Hall conductivity by the Kubo formula 
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Fig. 4. (Color online) Hall conductivity calculated at T = as 
a function of the electronic density n while assuming the four- 
sublattice all-out order [Fig. l(4a)] for (a) Jn = 0.22 and (b) 
Jll = 1.0. The shaded areas represent the regions where the 
four-sublatticc all-out order is realized at Jj^ = in Fig. 3(a). 

= E/ ^^d^nJ-^, (2) 

where fp is the Fermi distribution function, J^^'^ = 
(nQ,k| |?n/3k) [J^, is a current operator in the direc- 
tion V — x,y)^ and 'H jnck) — enk l^iak) with a de- 
generacy index a. Here, we take e'^/h = 1 (e is the 
elementary charge and h is the Planck constant). The 
second line is obtained by considering the non-Abelian 
connection £/nk = ipnkdtpnk, which is given by an M x 
M matrix-valued one-form associated with a multiplet 
■0nk = , InA/k)).^-^"^^) The integral is taken over 

the first Brillouin zone. For the four-sublattice all-out or- 
der, the multiplet with M = 2 describes the local SU (2) 
symmetry in the k-space owing to the band degeneracy 
associated with spins. When the Fermi level is in the 
gap, it is known that the Hall conductivity is quantized 
and given by the summation of Chern numbers for the 
occupied bands; cr^y = X]n(e k<eF) ^^ere the Chern 
number is defined as C„ — J^^^tr ds/nk/pisnkj/i'^Tri).^^^ 
Figure 4 shows the result of a^y obtained using eq. 
(2) while assuming the four-sublattice all-out order in 
the entire region of n for (a) Jh = 0.22 and (b) Jh = 
1.0. The shaded areas correspond to the regions where 
the four-sublattice all-out order is realized in the phase 
diagram at Jk = in Fig. 2(a). The results show that 
axy has a finite value in these chiral states. The Chern 
numbers associated with four bands are given by —1, 1, 
1, and —1 from the lowest band to the highest band, 
respectively. Thus, the Hall conductivity is quantized as 



o'xy = —e^/h at 1/4 filling and axy = +e^/h at 3/4 
filling when the system is a gapped insulator.^'*' This is 
observed in Fig. 4 at 3/4 filling for Jh = 0-22 [Fig. 4(a)] 
and at 1/4 filling for Jh = 1.0 [Fig. 4(b)]. 

To summarize, we have clarified that a four-sublattice 
noncoplanar magnetic order is stabilized around 1/4 fill- 
ing in the ferromagnetic Kondo lattice model on a trian- 
gular lattice. This phase is different from the 3/4- filling 
phase, which was predicted previously to appear owing 
to the perfect nesting of the Fermi surface. The region is 
wider than the 3/4- filling one, and includes both metallic 
and insulating states. It is further stabilized by the su- 
perexchange coupling between localized spins. The non- 
coplanar order retains a finite scalar spin chirality, re- 
sulting in the anomalous Hall effect; in particular, the 
Hall conductivity is quantized according to the Chern 
number in the insulating regions. Since our 1/4-filling 
phase is not directly related to nesting instability, it is 
expected to survive in a broad range of triangular-lattice 
systems, compared with the 3/4- filling phase. The sta- 
bility against a modification of band structure is a prob- 
lem for future study. It is also of interest to examine 
the stability against quantum and thermal fluctuations 
by, e.g., spin- wave analysis and Monte Carlo simulation. 
Our results will be relevant to spin-charge coupled phe- 
nomena in itinerant triangular-lattice compounds, such 
as delafossites. 
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